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Abstract. We report on an observation of the Anomalous X-ray Pulsar IE 1048.1-5937 performed with the XMM-Newton 
satellite. The phase averaged spectrum of IE 1048.1-5937 is well described by the sum of a power law with photon index 
~2.9 and a blackbody with temperature ~0.6 keV, without evidence for significant absorption or emission lines. The above 
spectral parameters do not vary during the phases corresponding to the broad pulse, while the off pulse emission shows a 
different spectrum characterized by a soft excess at energies below ~1.5 keV. The XMM-Newton observation and a re-analysis 
of archival BeppoSAX data, show that the spectral parameters and flux of IE 1048.1-5937 did not change significantly during 
observations spanning the last four years. All the data are consistent with a 2-10 keV luminosity varying in the range ~(5- 
7)x 10 33 erg s _1 (for a distance of 3 kpc). 
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1. Introduction 

The pulsar IE 1048.1-5937 is a persistent source of X-rays 
with flux of ^10~ n erg cm~ 2 s _1 and a significant modula- 
tion at a period of ^6.45 s (Seward, Charles & Smale 1986, 
Mereghetti 1995, Corbet & Mihara 1997, Oosterbroek et al. 
1998). It belongs to the small group of Anomalous X-ray 
Pulsars (AXP's, Mereghetti & Stella 1995, van Paradijs et al. 
1995). The nature of these pulsars, that show properties clearly 
different from those of the ^90 X-ray pulsars powered by ac- 
cretion from high mass companion stars, is currently unclear 
(see Mereghetti 2001, for a review). The lack of evidence for 
bright optical counterparts led to models for AXP's involving 
isolated neutron stars, but their typical values of P and P yield 
rotational energy losses inadequate to power the observed lumi- 
nosity. Although low mass companions cannot be completely 
ruled out (Mereghetti, Israel & Stella 1998), it has been pro- 
posed that the AXP's consist of isolated neutron stars accreting 
from residual disks (van Paradijs et al. 1995, Ghosh, Angelini 
& White 1997, Chatterjee et al. 2000). Another class of models 
is based on isolated, strongly magnetized (B>10 14 G) neutron 
stars, or "Magnetars" (Thompson & Duncan 1993, 1996; Heyl 
& Hernquist 1997), powered by the decay of their magnetic 
energy. 

Here we report on a new observation of IE 1048.1- 
5937 performed with the XMM-Newton satellite. The new data 
provide an accurate measurement of the flux, phase resolved 
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spectroscopy and a more precise localization of this pulsar. We 
have also reanalyzed the BeppoSAX data originally reported by 
Oosterbroek et al. (1998) in order to allow a more accurate 
comparison with the new observation. 

2. Observations and Data Analysis 

An observation of IE 1048.1-5937 was carried out with the 
XMM-Newton satellite (Jansen et al. 2001) on 2000 December 
28. Here we concentrate on the data obtained with the European 
Photon Imaging Camera (EPIC, Turner et al. 2001; Striider 
et al. 2001). The short duration of this observation (^8 ks), 
prevented a significant detection of IE 1048.1-5937 in the 
Reflection Grating Spectrometer (RGS). 

EPIC consists of two MOS CCD detectors and a PN CCD 
instrument. In total, the 3 EPICs provide > 2500 cm 2 of col- 
lecting area at 1.5 keV. The mirror system offers an on-axis 
FWHM angular resolution of 4-5" and a field of view (FOV) 
of 30' diameter. 

During this observation the medium filter was used. The 
MOS cameras were operated in Small Window mode to limit 
source photon pile-up effects and to improve the time resolu- 
tion. Due to its shorter frame integration time (73 ms), the PN 
could operate in the standard Full Frame mode. The net expo- 
sure time was 7135 s in each MOS camera and 4744 s in the 
PN. 

The EPIC data reduction was performed using version v5.0 
of the XMM-Newton Science Analysis System. An additional 
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Table 1. Summary of the Spectral Results^ 



Obs. 


Pow Law 


Absorption 


kT 


X 2 /d.o.f. 


2-10 keV Flux 


2-10 keV Flux 




photon index 


10 22 cm" 2 


(keV) 




absorbed (li) 


unabsorbed (erg cm~ 2 s -1 ) 


PN 


2.7±0.3 


0.9±0.1 


0.61±0.04 


132.1/129 


3.8 ±0.4 10" 12 


4.3 10~ 12 


MOS1 


3.1±0.4 


1.2±0.2 


0.63±0.07 


94.7/100 


4.2 ±0.5 10~ 12 


5.1 10~ 12 


MOS2 


3.3±0.6 


1.1±0.2 


0.71 ±0.07 


100.2/96 


4.0 ±0.5 10~ 12 


4.8 10" 12 


EPIC (6) 


2.9±0.2 


1.04±0.08 


0.63±0.03 


340.5/332 


4.1 ±0.4 10" 12 


4.9 10" 12 


BeppoSAX {c) 


3.3±0.4 


1.2±0.3 


0.62±0.04 


185.8/180 


5.9±0.5 10^ 12 


7.1 10" 12 



(a) Errors are at the 90% c.l. for a single interesting parameter 

(b) EPIC = PN+MOS1+MOS2 

Reanalysis of the observation reported by Oosterbroek et al. (1988) 
( d > The errors on the flux values are dominated by the systematic uncertainties (statistical errors are of the order a few percent) 




Fig. 1. EPIC spectrum of IE 1048. 1-5937 fitted with the sum of a power law and a blackbody (see the EPIC parameters in Table 
1). The upper spectrum is from the PN camera. 



data cleaning was required for MOS2 data in order to remove 
some residual low energy electronic noise. 

2.1. Spectral Analysis 

The target was clearly detected in the three EPIC cameras at 
coordinates R.A =10^ 50 m 06.3 s Dec. =-59° 53' 17" (J2000). 
Based on the current performance of the XMM-Newton satellite 
attitude reconstruction we conservatively estimate an uncer- 
tainty of ~4" on this position that is consistent with previous, 
less accurate measurements (Mereghetti, Caraveo & Bignami 
1992). 

The spectral analysis with the PN camera was based on 
counts extracted from a 100" x 100" box centered at the source 



position. This box was inside a single CCD chip. For the back- 
ground extraction we tried several source free regions close to 
the target, without finding significant differences in the results. 
All the spectra were extracted in the nominal 0.2-10 keV range 
and rebinned so that at least 20 counts were included in each 
energy bin and the channels did not oversample the PN energy 
resolution. 

We first analyzed spectra based on both single and dou- 
ble events. This gives a good statistics without significantly de- 
grading the spectral performances. 

Single component models, i.e. power law, thermal 
bremsstrahlung and blackbody, gave unacceptable results. A 
very good fit was obtained with the "canonical" AXP model 
based on the sum of a blackbody and a power law. The best fit 
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Fig. 2. Residuals of the PN single events spectrum 

parameters were a blackbody temperature kT^0.6 keV and a 
power law photon index a p h ~2.7. 

The two MOS spectra were based on all the counts of the 
small CCD windows (100"xl00") used in this observation. 
The background could only be extracted from the peripheral 
CCD chips, excluding small regions around sources and cor- 
recting the intensity for the vignetting effect that is significant 
in this region more than 5' offset from the telescope axis. We 
only considered the 0.2-10 keV energy range and applied the 
same rebinning as in the PN case. For both MOS spectra a com- 
bination of blackbody and power law gave the best fit result, 
with spectral parameters consistent with those obtained with 
the PN camera. 

We therefore proceeded to a joint spectral analysis of the 
three data sets, allowing only the relative normalization to vary. 
The results are reported in Table 1 for the blackbody plus power 
law model and the corresponding spectrum is shown in Fig. 1. 

The observed flux in the 2-10 keV energy range is 
~4.1 xl0~ 12 erg cm~ 2 s _1 in the two MOS cameras. The 
value obtained with the PN is about 10% smaller. Part of the 
discrepancy can be explained by the fact that in our analysis 
the PN "Out of Time" events are excluded. They are estimated 
to be the 6.2% of the total flux in Full Frame mode. The re- 
maining discrepancy is within the uncertainties in the relative 
calibration of the instruments. 

The best fit spectral parameters reported in Table 1 are con- 
sistent with those derived from previous ASCA observations 
(Paul et al., 2000), but they show some differences with re- 
spect to those observed with BeppoSAX by Oosterbroek et al. 
(1998), who found a harder powerlaw (a ph ^2.5) and a smaller 



absorption (N H ~ 0.45 x 10 22 cm -2 ). In order to check if this 
discrepancy is really due to a spectral variation we present a 
reanalysis of the BeppoSAX data in Section 2.4. 

2.2. Search for spectral features 

To fully exploit the energy resolution provided by the EPIC 
CCD detectors, we performed a spectral analysis by selecting 
only single events. Spectra based on single events, i.e. those for 
which the collected charge is not split between adjacent CCD 
pixels, provide the best energy resolution, at the expenses of a 
slightly reduced efficiency. 

To extract and rebin source and background spectra from 
the PN camera we followed the same procedure described in 
the previous section. The residuals of the best fit spectrum 
(power law plus blackbody) are shown in Fig. 2. Possible fea- 
tures are seen at energies of ~0.95 keV and in the region be- 
tween ^4 and 5 keV. 

Although the highest point in the low energy feature is at 
^3.5 o from the best fit model, the fact that its width is smaller 
than the instrumental resolution at this energy (FWHM ^60-70 
eV), suggests that it is most likely an artifact. 

The absorption at high energy has a width consistent with 
the instrumental resolution. It can be fit with a gaussian line 
centered at E=4.71±0.04 keV and with a=45 eV. 

None of these two possible features is confirmed by the 
MOS spectra (an excess is seen in the residuals at ^0.9 keV 
only in the MOS2 at -2 a). 
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Fig. 3. EPIC PN light curves of IE 1048. 1-5937 folded at the best period for different energy ranges 



2.3. Timing analysis and phase resolved 
spectroscopy 

The timing analysis was done by correcting the time of arrival 
of the counts to the Solar System barycenter and then folding 
them to several trial periods around the expected value. The 
best periods obtained for the PN and MOS data were respec- 
tively P=6.4526±0.0001 and P=6.4525±0.0001. 

The PN folded light curve is shown in Fig. 3 for different 
energy ranges. The pulse profile is characterized by a broad, 
nearly sinusoidal profile at all energies, as already noted in pre- 
vious observations with other satellites. The level of modula- 
tion was derived by fitting a constant plus a sinusoid to the 
normalized pulse profiles, yielding a pulsed fraction of —72% 
below 1.5 keV and ~94% at higher energies. 

To search for possible spectral variations as a function of 
the spin phase, we divided the data in four phase intervals cor- 
responding to the off-pulse (0.2-0.6), rising (0.6-0.8), top (0.8- 
1), and declining (0-0.2) parts of the light curve. The spectra 
of the four phase intervals were fitted with the power law plus 
blackbody model, allowing only the overall normalization to 
vary and keeping all the other parameters fixed at the best fit 
values of Table 1. The resulting residuals, see Fig. 4, clearly 
illustrate that the average spectral parameters provide a satis- 
factory fit to the three phase intervals covering the rise, peak 
and decay of the pulse. On the other hand, the spectrum dur- 
ing the phase interval 0.2-0.6 (corresponding to the minimum 
in the folded light curve) presents a significant excess below 
-1.5 keV. 



Despite the limited counting statistics, we attempted to 
model the spectral shape in the pulse minimum. A single power 
law provides an acceptable fit with parameters a p h~ 3.0 and 
N#~ 9x 10 21 cm~ 2 . Keeping the absorption fixed at the phase 
averaged value of 10 22 cm~ 2 , we obtain a p /j=3.2±0.1. 

2.4. Reanalysis of the BeppoSAX data 

We reanalyzed the BeppoSAX data of the May 1997 observa- 
tion already reported by Oosterbroek et al. (1998). The main 
difference with respect to the work of these authors is our es- 
timate of the background spectrum from a region of the field 
close to the target. Since IE 1048.1-5937 lies at low galac- 
tic latitude, where significant diffuse X-ray emission from the 
Carina nebula is present, the use of standard background files 
from BeppoSAX observations of extragalactic empty fields un- 
derestimates the background level (particularly at energies be- 
low —3 keV). 

Our extraction regions for the MECS (2-10 keV) were a cir- 
cle of radius 4' for the source and a circular corona from 280" 
to 560" for the background. The corresponding LECS (0.5-10 
keV) values were 8', 520" and 960". The best fit results to 
the joint LECS+MECS spectra for the power law plus black- 
body model are given in Table 1 . Both the absorption and the 
power law photon index are different from those reported by 
Oosterbroek et al. (1998) and consistent with the EPIC values. 
Thus we conclude that there is no evidence for a variation in 
the source spectrum between the two observations. The unab- 
sorbed flux corresponding to the best fit BeppoSAX parameters 
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Fig. 4. Results of the phase resolved spectroscopy. The figures show the ratios between the source spectra in the indicated phase 
intervals and the average spectrum. 



is 7.1 x 1CT 12 erg crrr 2 s _1 (2-10 keV), about 45% higher than 
that observed with XMM-Newton. 

3. Discussion 

The XMM-Newton observation of IE 1048.1-5937 provides 
the first clear evidence for a spectral variation as a function 
of the pulse phase in this source. The spectrum shows a "soft 
excess" (with respect to the phase averaged spectrum) during 
the pulse minimum, while no significant variations occur along 
the pulse rise and decay. This behavior is different from that of 
other AXP's which show energy dependent pulse shapes like 
1RXS J 170849-4009 (Israel et al. 2001, Gavriil & Kaspi 2001) 
and 4U 0142+61 (Israel et al. 1999). 

Although the spectra of all AXP's are well described by the 
sum of a power law and a blackbody, there are no compelling 
reasons to attribute this spectral shape to the presence of two 
physically distinct emitting components (some suggestions in 
this sense were made, e.g., by Ghosh, Angelini & White 1997). 
In fact, as noted by Ozel et al. (2001), the contribution of the 
blackbody component to the total flux is a strong function of 
the energy, while the pulsed fraction at different energies does 
not change significantly. In the case of two physically distinct 
spectral components this would require an ad hoc coupling in 
the pulse profiles of the blackbody and power law components. 

The fact that the pulsed emission from IE 1048. 1-5937 can 
be fit by the same spectral shape at all the phases also sug- 
gests that it originates from a single physical component, with 
a spectrum more complex than a pure blackbody. Theoretical 



AXP's spectra have been recently computed by several authors 
in the context of models of thermal emission from strongly 
magnetized neutron stars. These models predict modifications 
to pure blackbody-like spectra, in the form of hard power-law 
tails (Ozel 2001, Perna et al. 2001) and/or narrow spectral fea- 
tures (Zane et al. 2001). Further observations with better statis- 
tics are clearly required to confirm the possible features in the 
spectrum of IE 1048.1-5937 reported in section 2.2. 

The period value found in the XMM-Newton observation 
is consistent with the overall spin-down trend at an average 
P~1.6xl0- n s s" 1 that IE 1048.1-5937 has maintained in 
the last few years. Large variations in the average spin-down 
rate were detected in the past. In particular, between June 
1992 and March 1994 (Mereghetti 1995, Corbet & Mihara 
1997), P~3.3xl0~ n s s" 1 was about a factor 2 higher than 
that measured earlier and after 1996. More recent observations 
have clearly shown that the level of rotational irregularities in 
IE 1048.1-5937 is higher than in other sources of this class 
(Paul et al. 2000, Baykal et al.2000, 2001; Kaspi et al.1999, 
2001; Gavriil & Kaspi 2001). 

In principle, the study of the spin-down evolution in 
AXP's could discriminate between accretion-based and mag- 
netar models. If AXP's are powered by accretion (either from 
a companion star or from a residual disk) one would expect 
significant P fluctuations superimposed on the long term spin- 
down trend, possibly correlated with luminosity variations, due 
to irregularities in the accretion flow. On the other hand, the pe- 
riod evolution of a magnetar should be much more regular, with 
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Fig. 5. Flux and spin period evolution of IE 1048.1-5937 during the years 1996-2000. The top panel reports the 2-10 keV flux 
corrected for absorption as measured with BeppoSAX, ASCA and XMM-Newton. The different symbols in the bottom panel are 
as follows. Triangles: BeppoSAX and XMM-Newton (this work); Squares: RossiXTE (Mereghetti et al. 1998, Baykal et al. 2000); 
Diamonds: ASCA (Paul et al. 2000); Solid lines: RossiXTE (Kaspi et al. 2001). 



the possible exception of "glitches" as observed in young radio 
pulsars. 

It is unclear whether the spin-down variations seen in 
IE 1048. 1-5937 are related to changes in its X-ray flux, partic- 
ularly for observations carried out before 1995. Oosterbroek et 
al. (1998) compared the flux of IE 1048. 1-5937 measured with 
BeppoSAX to all the values obtained in the previous years with 
the Einstein, EXOSAT, ROSAT and ASCA satellites. Although 
the flux values span almost an order of magnitude, the system- 
atic uncertainties are difficult to evaluate, since some of these 
measurements referred to different energy ranges. Furthermore, 
the data from non-imaging detectors might be contaminated by 
the presence of the strongly variable source Eta Carinae within 
the field of view. 

By means of a systematic program aimed at phase-coherent 
timing of IE 1048.1-5937 with the RossiXTE satellite, Kaspi 
et al. (1999) showed the presence of significant P changes 
between 1996 and 2000. Their ephemerids are shown by the 
solid curves in Fig. 5. Unfortunately, due to difficulties in the 
background subtraction, the RossiXTE non-imaging observa- 
tions can measure only the intensity of the pulsed flux, which 
was found consistent with a constant value. Our XMM-Newton 
and BeppoSAX results, compared with the ASCA ones (Paul et 
al. 2000), show that all the measurements obtained since 1994 
with imaging detectors are consistent with an unabsorbed 2-10 
keV flux within the range ~(5-7) x 10~ 12 erg cm~ 2 s _1 . 



Our period measurement is inconsistent with an extrapo- 
lation of the phase-coherent solution reported by Kaspi et al. 
(2001). This is a further demonstration of the high level of tim- 
ing noise present in IE 1048.1-5937 . 

4. Conclusions 

The large collecting area of the XMM-Newton telescope has 
provided the first significant evidence for spectral variations as 
a function of the spin phase in IE 1048.1-5937 . Contrary to 
other AXP's, the spectral shape changes only in the off-pulse 
emission, while no significant variations are seen during the 
broad pulse. 

No convincing evidence was found for narrow spectral fea- 
tures. Longer XMM-Newton observations (also exploiting the 
RGS instrument) and improved calibrations of the instruments 
are required to verify (or not) the structures marginally detected 
in the EPIC PN camera. 

The XMM-Newton observation, as well as our re-analysis of 
the archival BeppoSAX data, show that the spectral parameters 
and flux of IE 1048. 1-5937 did not change significantly during 
observations spanning the last four years. The fact that all the 
observations of good quality obtained with imaging detectors 
and covering energies up to ^ 10 keV are consistent with a 2-10 
keV luminosity varying in the range ~(5-7) x 10 33 erg s _1 (for 
d=3 kpc), casts some doubts on previous reports of variability 
up to a factor 10. It is unclear whether the large timing noise 
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of IE 1048.1-5937 can be ascribed to variations in the torque 
exerted on the neutron star by such relatively small fluctuations 
in the mass accretion rate. 

Finally, we have provided a reduction of a factor ~15 in the 
area of the error box of IE 1048. 1-5937 . This will undoubtedly 
help in the search for an optical counterpart of this intriguing 
X-ray source. 
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